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In the junctions with a strong Josephson coupling and a pronounced interfacial pair breaking, the 
magnetic interference patterns of the Josephson current are shown to differ substantially from the 
standard Fraunhofer shape. The Fraunhofer pattern occurs, when Josephson couplings are weak. 
The narrow peak of the critical current, centered at the zero magnetic field, and the suppressed hills 
at finite field values are the characteristic features of the non- Fraunhofer magnetic field modulation 
of the critical current, obtained in this paper. 
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The Josephson critical current, as a function of the 
applied magnetic field, manifests a characteristic modu- 
lation, which originates from the fundamental quantum 
coherence of the superconducting leads. Under standard 
conditions the modulation is described by the Fraunhofer 
diffraction pattern, which has represented a benchmark 
for the magnetic field modulation of the critical current 
since the early days of the discovery of the Josephson 
effect. [H-Q Distortions of the Fraunhofer shape in tun- 
nel junctions stem from inhomogeneities or from vari- 
ations of geometry of the fabricated interfaces. 0, H[ 
In contemporary studies of the magnetic field induced 
interference patterns, single junctions are usually fabri- 
cated as a sequence of different regions, e.g., of and 
7r regions. In such junctions the patterns can substan- 
tially deviate from the standard shape, due to their phase 
sensitivity. The corresponding measurements constitute 
the basis of the Josephson intcrfcromctry method [Ml, 
which has been successfully used to identify symmetries 
of, e.g., the d-wave order parameters of cuprate super- 
conductors [5k7 |, the complex p-wave superconductivity 
in Sr2RuC>4 [TJ, 1 1 311 and the Ei u pairing of the super- 
conducting UPt3 [14j . Single Josephson junctions with 
inhomogeneous magnetic intcrlayers, fabricated as a se- 
quence of and ir regions, also exhibit the interference 
patterns, which qualitatively differ from the Fraunhofer 
ones. [liOl 

This paper theoretically studies single junctions with 
homogeneous thin interfaces and focuses on the effects of 
the anharmonic phase dependence of the Josephson cur- 
rent on the magnetic field induced modulation. The mod- 
ulation of the Josephson current will be described within 
the Ginzburg-Landau (GL) theory. The effects will be 
shown to become pronounced in the junctions with a 
strong Josephson coupling, where the modulation differs 
heavily from the Fraunhofer diffraction pattern. The re- 
sults obtained could be, in particular, implemented in 
the junctions involving unconventional superconductors 
and/or magnetic interfaces. An intensive interfacial pair 
breaking is crucial for the systems in question. It sub- 
stantially suppresses the critical current, therefore main- 



taining the planar junctions with a pronounced Joseph- 
son coupling as the weak links with strongly anharmonic 
current-phase relations. [l9| 

Assume the usual form of the GL free energy, which 
applies, for example, to s-wave and c? a ,2_ l( 2-wave junc- 
tions. Let the symmetric junctions have a spatially con- 
stant width much less than the Josephson penetration 
length. The superconducting electrodes Si and S r are 
supposed to be thick compared to the magnetic pene- 
tration depth A, while a homogeneous plane rectangu- 
lar interlayer at x = is assumed to be of zero length 
within the GL approach (see Fig.[l|. The interface con- 
tribution to the GL free energy incorporates both the 
Josephson coupling of the two superconducting banks 
g.j\^ + - *-| 2 and the term g{\^+\ 2 + |*-| 2 ), which is 
responsible for the interfacial pair breaking or pair pro- 
ducing. It will be convenient to use the dimensionlcss 
Josephson coupling constant gi = gj£(T)/K, the inter- 
face pair breaking parameter g$ = g£(T)/K and their 
combination gb(x) — (gs + 2<7fsin 2 -|), where £(T) is 
the temperature dependent superconducting coherence 
length, K comes from the bulk gradient term i^lV^I 2 
and x is the phase difference of the order parameters 
across the interface. 




FIG. 1. Schematic diagram of the junction. 

Microscopic estimations show that the range of vari- 
ations of the parameters gi and gs (and, therefore, of 
<7b) is quite wide and, in particular, includes large values. 
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Til 20j The limit 5,5 — > +00 corresponds to the order pa- 
rameter vanishing at the boundary. In highly transparent 
junctions ((1 -D) < 1) the parameter gi oc (1 — I?) -1 
can be arbitrarily large. Here I? is the transmission co- 
efficient. For small and moderate transparencies in dirty 
s-wave junctions gi can vary from vanishingly small val- 
ues in the tunneling limit to those well exceeding 100 
near T c . For the pair breaking surfaces (<7& > 0) arbi- 
trarily large values of gb are admissible in the GL theory, 
since they correspond to spatial variations of the order 
parameter on a scale > £(T). This is in contrast to the 
pair producing interfaces gb < 0, where the supercon- 
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ductivity is locally enhanced (see, e.g., 
the boundary a strong enhancement can induce the scale 
substantially less than the coherence length of the leads. 

The self-consistent description of the Josephson cur- 
rent in planar junctions, recently developed within the 
GL theory in Ref.[l9|, takes into account the pair break- 
ing effects produced by the phase-dependent Josephson 
coupling as well as by the current and by the interface 
itself. In particular, the following anharmonic current- 
phase relation has been obtained in the absence of the 
applied magnetic field: 
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Here jd P is the depairing current in the bulk. 

Equation ([I]) describes the current behavior almost 
perfectly, if j < 0.7jd P . This concerns, in particular, 
the current at ge < 1 for any gs, or at gs > 1 for any 
gi. For j > 0.7 jdp equation ([T]) reasonably approximates 
the exact numerical solution, with the deviations not ex- 
ceeding 10%. [lj| The solution ([T]) applies only to the 
pair breaking interfaces, i.e., to gb(x) ^ 0- Otherwise, 
the parameters in equation ([T]) can take arbitrary values. 
In order to describe the current flowing through the pair 
producing interface, one should change the sign before 
the square root in equation ([I]). Depending on gn and 
gs, such current could exceed jdp and, hence, destroy the 
superconductivity in the bulk of the leads. From now on 
only the pair breaking interfaces with gi,gs > will be 
considered. 

For a pronounced interfacial pair breaking g 2 3> 1 ex- 
pression ([T]) for the supercurrent is simplified and reduced 
to 
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The corresponding critical current j c = [3V3ge/4gs(gs + 
25^)]jdp *C jdp is always small, at arbitrary ge. The va- 
lidity of the condition j <C jdp is the evidence of a weak 
link. Strongly anharmonic current-phase relation shows 



up in ([2j for g 2 ^> g 2 ^> 1, while in the case gt <C gs the 
first harmonic dominates the current. 

Now let the magnetic field be applied to the junction 
along the z axis: B(x) = B{x)e z (see Fig. [I]). It is as- 
sumed that the field is not too strong, substantially less 
than the critical fields of the leads. The well known gauge 
invariant result states that the supercurrent density in 
the presence of the field is described by its zero-field ex- 
pression j(x), if the spatially dependent phase difference 
across the interface x(y) = X + / L y )(& /&o) , influ- 
enced by the magnetic flux <fr through the junction, is 
substituted for x- S HI Here $0 = 7rftc/|e| is the su- 
perconducting flux quantum. Considering the total su- 
percurrent through the junction is of interest, the spa- 
tially modulated supercurrent density should be averaged 

over the interface area: / = -r-J 2 Ly j[x+ j^$r ) dy = 

y \ y J 

t|^$ fp + j(ip)d(p. Here a rectangular plane interface is 
supposed to occupy the space (—L y /2, L y /2) along the 
y axis. The integration limits (p± arc defined as <p± = 
X±tt($/$ ). 

A macroscopic scale of the modulation period Ly = 
ir£ 2 3 /2\ = $ /2i?(0)A allows one to consider, for a 
given y, local densities per unit interface area of ther- 
modynamic potential and of the supercurrent, which 
satisfy the relation j(x) = 2 j^-- d ^f2o(x)- Here £3 = 
(Hc/leBfi)}) 1 / 2 is the magnetic length. With this rela- 
tion, the Josephson current averaged over the interface is 

M 



(3) 



Thus, the phase dependent thermodynamic potential of 
the junction in the absence of the field also describes the 
magnetic interference pattern in the junction. 

Equation §3§ will be now used to establish the modula- 
tion of the anharmonic Josephson current by the applied 
magnetic field. Interestingly, for a strong interfacial pair 
breaking g 2 3> 1, when the supercurrent in the absence 
of the field is described by © , the averaged anharmonic 
current allows not only a numerical study but also an 
analytical description. Indeed, integrating @ over the 
phase difference and using ([3]), one obtains the magnetic 
field dependent current-phase relation for the averaged 
anharmonic supercurrent: 
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(a) tunnel junctions (gi <C 1) (b) gi = 10 (c) gi = 100 

FIG. 2. The current-phase relations in junctions with the pair breaking parameter gs = 4, taken at several values of the 
magnetic flux $: (1) $ = (2) $ = 0.3$o (3) $ = 0.6$o (4) $ = 0.9$o . Here <E>o = nhc/\e\ is the flux quantum. Different 
panels correspond to different Josephson coupling constants. For each coupling constant the current is normalized to its critical 
value in the absence of the field. 



The corresponding critical current is 
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where A = 2 (gs + gi) t 

Equations (j4|) and (0, as well as the plots shown in 
FigS-HJ [3] and 0] represent the central results of this pa- 
per. Fig. [2] displays the current-phase relations obtained 
numerically with the exact self-consistent calculations 
within the GL theory for = 4 and for various values 
of gi and 5>. The corresponding equations used for this 
purpose will be outlined in the end of the paper. The ap- 
proximate results based either on equations ([T]) and (|3|), 
or on the analytical formula (|4|, all lead to the curves, 
which are in excellent agreement with and practically in- 
distinguishable from the exact ones shown in Fig. [2 Sim- 
ilarly, the analytical expression agrees well with the 
numerical results for the magnetic interference patterns 
shown in Fig.[3j 

The modulations obtained substantially differ from 
the Fraunhofcr pattern in the junctions with the 
strongly anharmonic current-phase dependence, which 
takes place for large Josephson couplings gf 3> fff 3> 



1. Figs. 2(b) 2(c) demonstrate that in such junctions 



the magnetic field substantially modifies the anharmonic 
structure of the current-phase relations. This is in con- 



trast to tunnel junctions (see Fig. 2(a)), where the sinu- 
soidal shape of the phase dependence is not distorted by 
the magnetic field. 

An important feature of the non-Fraunhofcr patterns, 



obtained under the conditions gi 3> gs, g% 3> 1, is 
that the zero-field peak gets narrower with increasing 
gt,. As seen in Fig. 2] and as follows from equation 
the half width at the half of the peak is small for large 
values of gi/gs and can be approximated as (A$/$o) ~ 




FIG. 3. The normalized critical current 7 C / I c o as a function of 
the magnetic flux, which pierces the junctions with different 
Josephson coupling constants: (1) small Josephson couplings 
gi <C 1) described by the Fraunhofer pattern (2) ge — 10 
(3) gi = 30 (4) ge = 100. The interfacial pair breaking 
parameter is gs — 4; 7 c o is the zero-field critical current for 
the corresponding junctions. The curves have been plotted 
based on equations |T]) and 
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FIG. 4. The half width at the half of the zero-field peak as 
a function of the Josephson coupling constant ge, according 
to equation ([5]). The curves are plotted for various values of 
the interfacial pair breaking parameter gs: (1) gs = 4 (2) 
gs = 20 (3) gs = 50 (4) g s = 100 (5) gs = 500. 



l.35gs/gt -C 1. In the opposite case gi/gs -C 1 the quan- 
tity (A$/$o) is close to its Fraunhofer constant value 
w 0.6 and weakly depends on gg. Also, the hills at the 
finite fields shown in Fig. [3] are significantly suppressed 
as compared to the Fraunhofer ones, in measure of the 
same small parameter gs/ge "C 1 and up to a logarithmic 
factor. The larger the Josephson coupling, the stronger 
the suppression of the critical current by the magnetic 
field within each interval n$o < < (n + l)$o (n = 
0, ±1, ±2, . . . ). The Fraunhofer pattern is represented 
by the curve 1 in Fig. [3] It appears for ge <C gs, when 
the first harmonic dominates the current. Then one gets 

from (0 I c ($) = I c0 |sin (7r<f>/(f>o) | / (7r|$|/<f> ) . where 

IcO = JdpSVSge/ig 2 is the corresponding critical current 
in the absence of the magnetic field. 

Let us return now to a brief description of the approach 
used for plotting Fig. [2] and based on the results of Ref. 
[l9| . For the order parameter f(x)e lx ^ normalized to 
/ = 1 in the bulk without supcrflow, one writes the first 
integral of the GL equation in the presence of the super- 



current [27 1 
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Here x = aj/f(T), j = j/j dp = - (3^3/2) (d X /dx)f 2 and 
foo is the asymptotic value of / in the depth of the leads. 

The boundary conditions for / and the expression for 
the supercurrent in the symmetric junctions are 



±5b(x)/o, 3 
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One puts x — in and, using (|7j). eliminates the 
current and the first derivative of the order parameter. 
This results in a biquadratic relation between /g and 



f^. Then, using the asymptotic formulas in the bulk 
j = (3V3/2)v s (l-v 2 s ),f a 



1 



and equating the 
current to that in {7J with /q = (1 — v 2 )a, one obtains 
v s = age sin \- As v s , fa, foo and j are now expressed 
via the only variable a, one gets from the biquadratic 
relation the fourth-order polynomial equation for a [l^l 

2g 2 ( X )a - (1 - a) 2 [l - a(a + 2)g 2 sin 2 X ] = 0. (8) 

Eq. (O is exact within the conventional GL ap- 
proach with the boundary conditions ([?])■ Calculating 
a{9t,9s,X), one S ets j(x) and finds numerically J7 (x)- 
Then, from equation with the calculated i?o(x), one 
obtains the self-consistent current-phase relations shown 
in Fig.CU 

In conclusion, it is demonstrated that the anharmonic 
phase dependence of the Josephson current, taking place 
in the junctions with a strong Josephson coupling and 
a pronounced interfacial pair breaking, has a profound 
influence on the phase sensitive magnetic interference 
patterns. The distinctive features of the non-Fraunhofer 
patterns uncovered, are the narrow peak of the critical 
current, centered at the zero field, and the suppressed 
hills at finite field values. The results obtained could be 
implemented in the junctions involving unconventional 
superconductors and/or magnetic interfaces. 
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